Accelerating Universes from Short-Range Interactions. 
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We show that short-range interactions between the fundamental particles in the universe can 
drive a period of accelerated expansion. This description fits the early universe. In the present day 
universe, if one postulates short-range interactions or a sort of "shielded gravity", the picture may 
repeat. 
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There is no doubt that the inflationary paradigm has 
played a central role during the last several decades in our 
understanding of the early universe QJ. While most of 
the authors rarely question the initial acceleration itself, 
the underlying cause of the acceleration remains rather 
obscure. It is often assumed that the initial accelerated 
expansion is driven by a sort of a self interacting scalar 
field , the so-called inflaton. Much work has been done 
during the last twenty years studying the possible form 
of the potential for the scalar field which best fits within 
the cosmological model Q. Recently, other proposals 
have been put forward where the inflaton takes the form 
of a non-canonical scalar field, also known as K-field Q. 
Needless to say that the inflaton is rather an effective field 
which most probably does not correspond to any funda- 
mental particle, but could be looked at as an effective 
description of an underlying physical theory jjj. 

A related issue is the state of the matter at the very 
high densities which prevail in the early universe. Pre- 
sumably, the fundamental theories such as Superstrings, 
M-theory, etc. ;(j might guide one as to what symmetries 
and laws to expect when the densities, velocities and en- 
ergies approach those present near the Big Bang. Certain 
clues about the behavior of the ultradense matter can be 
also obtained in the framework of Quark-Gluon Plasma 
theory pj, which could be handy to model the description 
of the primordial matter in the early universe. This area 
of research is now under intensive theoretical and exper- 
imental study, mainly at RHIC and at CERN, with the 
surprising new observation that in this extreme state the 
matter seems to behave almost as a non-viscous perfect 
fluid 0|. At this stage, however, there is little one can 
say with a certain degree of rigor about the properties of 
matter nearly the initial singularity. 

This brings us to the following question: Is there a way 
to obtain inflationary solutions in a highly dense universe 
within a so-to-say "conventional physics"? Physics, which 
one might label as fairly realistic and, to some extent, 
generic. The answer to this question is the main purpose 
of this Letter, and follows to be affirmative. 

We will be considering a universe filled with an ultra- 
dense matter and modeling the interaction between the 
particles by a short-range attractive force. For compu- 



tational purposes we will be using the Yukawa type po- 
tentials, yet, our discussion is generic and applies to any 
attractive short-range interaction, as is explained below. 
Moreover, our findings hold through as well in the case 
of the diluted matter, but in this case, one would need to 
postulate some unknown short-range interaction, or as- 
sume that the Newtonian gravity in the expanding uni- 
verse is somehow "shielded". We will comment on these 
issues before closing. Our main result is, that describing 
the matter as a fluid of interacting particles via short- 
range attractive forces, the cosmological models undergo 
a phase of accelerated expansion. 

The introduction of short-range forces between the 
particles leads to the following equation of state for the 
matter: 



p = mon — An , p 
or equivalently, in a p(p)-form, 
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Here A is a term depending on the interaction, mo is 
the rest mass of the fundamental particle and n is the 
particle number density. The conventions and units we 
use are c = 8irG = h = 1 and the metric has a signature 
(-,+,+,+)■ 

This, somewhat unusual in cosmology equation of 
state, is the one which could drive the accelerated ex- 
pansion in the early universe. We see that if one neglects 
the interaction between particles (A = 0), or consider 
a very dilute matter (n — ► 0), the equation of state (TIJ 
reduces to the standard dust equation p — m^n, p = 0. 
Here we impose also the positivity of the energy density 
(n < mo/A). 

To obtain the above equation of state (see for exam- 
ple Q), we consider a system of N identical interacting 
particles placed at points 1, (i = 1, 2, ...,N). The inter- 
action between the particles is modeled by a potential 
V (hj), where we have defined kj = |1$ — L,-|, We further 
neglect the temperature effects by assuming that the par- 
ticle masses and the strength of the interaction are much 
larger than the effects of the temperature, and write the 
energy of the system as 

N 
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Dividing by the volume and assuming that the system 
is homogeneous and isotropic, at least at the scale of 
the interaction, one readily obtains the energy density 
as a function of n as given in the equation (QJ. The 
pressure, then, is computed using the thermodynamic 
relation p — n {dp/dn) — p. The interaction term A can 
be evaluated passing as usual to the continuous limit and 
replacing the sum by the integral: 
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V(l)l 2 dl. 



(2) 



Note, that the sign of the parameter A defines the sign 
of the pressure, and therefore, in order to obtain an ac- 
celerated expansion, the interaction should be attractive. 

A typical way of modeling a short range interaction 
is via the Yukawa potential V(l) = —g 2 e~ t * l /l. Here g 
is the coupling constant of the theory and p the mass 
of the boson mediating the force, whose Compton length 
lo = 1/p defines the range of the interaction. The integral 
of the equation (J2J, therefore, is unaffected by the upper 
limit at infinity. Performing the integral, we obtain: 



A = 



2irg 2 
V 2 



It is important to say that the convergence of the above 
integral basically defines as to whether the interaction is 
short- or long-range, and whether the system of parti- 
cles may be treated by conventional thermo- and hydro- 
dynamics. If the potential energy due to the interaction 
between the particles falls faster than r -3 , the force is 
short-range, the above integral converges and the usual 
thermodynamics apply. Such is the case for the Yukawa 
potential or, for example, the Van der Waals forces, etc. 
If, however, the potential energy does not fall as fast as 
r~ 3 , as for example in the case of Newtonian gravity one 
deals with long-range forces, and the dynamics of the sys- 
tems governed by those may be quite involved [TQl . There 
is little doubt, however, that the forces between the fun- 
damental particles in the early universe are short-range, 
therefore providing a broad motivation for our further 
analysis. 

We have found it useful to express the pressure, using 
the thermodynamic relations, in terms of the enthalpy 
per particle h: 

P(h) = -^(h-m ) 2 , (3) 
where the enthalpy is given by 

h=^- = m - 2An. (4) 
on 

We further assume that the ultradense matter is de- 
scribed by an isentropic irrotational perfect fluid, and to 
proceed, we give the action from which the equations of 
the fluid motion are derived fill ll^ : 



s= f dW=s{p(|v|)- (j£ 
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Here the current V 1 , usually known as the Taub current, 
is given by — hu^ (|V| = h), where u^, the 4- velocity 
of the fluid, verifies u^u^ = — 1. The dynamical variables 
are g^ v , and 4>, and the following equations of motion 
result: 



-h- l 6. 



with the stress-energy tensor given by 
2 6S 



rpfj.il 



dp 
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The usual thermodynamic relation between the density 
and the pressure, p = nh — p (with n = (dp/dh)), al- 
lows one to cast the energy-momentum tensor into the 
standard perfect fluid form: 

T» v = (p+p)u»u v +pg^. 

The equation © is the Clebsch decomposition of the 4- 
velocity for an irrotational flow, whereas the equation © 
is the conservation law for the particle number. Equation 
(HJ shows that the scalar field <j> plays the role of a velocity 
potential. 

The identity u^u^ = — 1 and the equation © lead to 
the following expression for the enthalpy in terms of the 
derivatives of the velocity potential: 



(8) 



Taking into account the Clebsch decomposition of the 
4- velocity 101, and introducing the new variable X = 
— \g^ v §^4>^ = h 2 /2, we obtain the on-shell expression 
for the action: 



Jon — shell 



-gF{x), 



where in the final expression we have defined 



p(h) = ph 



F(X). 



(9) 



(10) 



Expression © gives the action for an irrotational per- 
fect fluid. This functional takes the form of a non- 
canonical scalar field action and is analogous to the one 
often used in K-essence cosmology (known as purely ki- 
netic K-essence) (j,ll2,ll3illJ|- The Lagrangian density 
is obtained from the equation of state which relates the 
pressure with the enthalpy 111 Oil , and the scalar field is 
the velocity potential of the fluid, which plays the role 
of the inflaton field. For completeness, the density and 
particle number in terms of the variable X are given: 

p = 2XF'(X) — F(X), n = ±V2XF'(X), (11) 

where F'(X) denotes the derivative of the function with 
respect to its variable. 
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Figure 1: The figure shows the shape of the function F in 
terms of the kinetic term X. The upper branch, above the 
point q, represents the standard matter, whereas the lower 
branch dives into the phantom domain. The solid line repre- 
sents the values of X for which inflation occurs. The values 
of the different points are: X s — X p = b 2 /4a 2 , X r = b 2 /16a 2 
and F(0) = -b 2 /4a. 



The above formalism together with the equation of 
state in the form 10 and the assumption that one deals 
with an irrotational perfect fluid, gives the following mat- 
ter Lagrangian: 



F±(X) = -aX±b4x , 

4a 

where the parameters a and b are given by: 
1 m 



(12) 



2A' 



V2A' 



The ± sign in the equatio ns a ppears due to the definition 
of the enthalpy (h = ±\/2X) and separates the matter 
into an ordinary one, in the case of positive enthalpy, 
and the so-called phantom matter (p + p < 0) in 
the case where the enthalpy is negative. It is interesting 
that in the case of the perfect fluid, the phantom and 
the ordinary matter have a simple physical separation 
depending whether the enthalpy per particle is negative 
or positive respectively. 

The behavior of the function F(X) is depicted in the 
figure 1. Note, that the function F(X) is defined by the 
physical parameters of the system: mn and A and is not 
assumed ad hock as is usually done in the context of the 
K-essence cosmology. The positivity of the energy and 
of the number of particles imposes X < X s for both the 
ordinary and the phantom branch. The condition of the 
accelerated expansion p + 3p < is always verified in the 
phantom case, whereas in the case of the ordinary matter 
one needs X < X r for inflation. 

We now consider an isotropic and homogeneous uni- 
verse, modeled for simplicity by the spatially flat 



Friedmann- Robertson- Walker metric: 

ds 2 = -dt 2 + R 2 {t) [dr 2 + r 2 (d9 2 + sin 2 6dip 2 )] . (13) 

R{t) is the scale factor and (t, r, 9,ip) the comoving co- 
ordinates. The physical distances in such a universe are 
given by I = R(t)r. The dynamical equation for the field 
X can be obtained either from the particle conservation 
equation (jSJ or by the direct variation of the on-shell 
action and results in: 



[F'(X) + 2XF"(X)] X + 6HF'(X)X = 0. 



(14) 



Here, H = R/R is the Hubble parameter. Solving this 
expression for X, we obtain: 



X 



6HF'(X)X 
' F'{X) + 2XF"{X)' 



(15) 



The equation l(T5)l basically determines the behavior of 
the cosmological model We assume we deal with 

an expanding universe by fixing the sign of the Hubble 
parameter to be positive. Since the field X is positive 
by definition, the sign of X is determined by the signs 
of F'(X) and F'{X) + 2XF"(X). We can check that 
F'(X) + 2XF"{X) = -a for both branches, whereas 
F'(X) is positive on the positive branch of our picture, 
and negative on the negative one. This means, that what- 
ever the initial value of X is, one always finishes with the 
dust behavior, i.e. the pressure vanishes. If the initial 
value of the field X falls into the phantom branch, X 
decreases to zero, passes to the standard matter branch 
and then increases, till finally the model becomes pure 
dust near X s . The interesting point is, however, that 
the short-range interaction between the particles leads 
to a period of accelerated expansion naturally finished at 
X r , which is an exit point towards the dust-like universe 
(grateful exit). 

The picture that emerges, is therefore, as follows. The 
short-range attractive forces in the early universe intro- 
duce a measure of negative pressure (tension) and could 
be responsible for the early universe inflation. The in- 
flation stops naturally at the exit point defined by the 
parameters of the theory. 

One should not expect from this simple model, as it 
stands, a reasonable amount of inflation. Indeed, our es- 
timates do not seem to produce anything close to 60-70 
e-foldings one would like to have for an effective infla- 
tion. The point is that the model is too simple to aco- 
modate the interplay between the scales of the particle 
horizon and the interaction range. Moreover, in the case 
of the very early universe there will be instants where 
the horizon size would be even smaller than the range 
of the interaction. Then, the interaction term A will be 
time dependent, due to the fact that the upper limit of 
the integral in the equation (01 will be bounded by the 
size of the time-dependent particle horizon. This, intro- 
duces interesting physics and enhances the inflation, on 
one hand, but complicates the equations of the dynami- 
cal evolution on the other. We hope to be able to report 
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the results of the numerical studies of these equations in 
the near future. 

It would be extremely interesting if a similar mecha- 
nism could explain as well the present day acceleration 
[lfj . We see here two possible alternatives. First, would 
be to postulate the existence of fundamental particles 
which dominate the universe and interact via short-range 
forces m the sense that the expression converges. 
We will not speculate about this possibility here, how- 
ever, it should be pointed out that interesting models 
have been recently proposed in which the late time ac- 
celeration of the universe is obtained with the use of a 
Van der Waals equation of state ^3], and, as is known 
from statistical mechanics, this kind of equation of state 
appears when one takes into account the interaction be- 
tween the particles. 

Another possibility, which looks more appealing to us 
and is connected again to the interplay between the hori- 
zon size and the range of the interaction, is to consider 
the Newtonian gravity between the galaxies, or clusters, 
as the dominant contribution to the matter density and 
the pressure in the present day universe. One can think 
of the Newtonian gravity as the interaction between the 
fundamental particles in the universe. "Averaging" on a 



scale comparable with the cluster scale one can consider 
this term as an effective part of the energy-momentum 
tensor which drives the expansion. Now, the Newtonian 
gravity is a long-range force. Nevertheless, in an expand- 
ing past singular universe, one may apply a natural cut 
off to the integral (j2j to evaluate the interaction term A. 
Since no interaction (including gravity) acts beyond the 
particle horizon, the latter may serve as the interaction 
range scale for the gravity. Effectively, then, the grav- 
ity would become shielded by the horizon in an infinite 
past-singular expanding universe and could produce the 
necessary negative pressure to accelerate the universe. 
We leave, however, this and the related problem of the 
time-dependent interaction term for future report. 
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